Steatotic liver disease induced by TCPOBOP-activated hepatic constitutive androstane receptor: primary and secondary gene responses with links to disease progression.

Constitutive Androstane Receptor (CAR, Nr1i3), a liver nuclear receptor and xenobiotic sensor, induces drug, steroid and lipid metabolizing enzymes, stimulates liver hypertrophy and hyperplasia, and ultimately, hepatocellular carcinogenesis. The mechanisms linking early CAR responses to later disease development are poorly understood. Here we show that exposure of CD-1 mice to TCPOBOP, a halogenated xenochemical and selective CAR agonist ligand, induces pericentral steatosis marked by hepatic accumulation of cholesterol and neutral lipid, and elevated circulating alanine aminotransferase, indicating hepatocyte damage. TCPOBOP-induced steatosis was weaker in the pericentral region but stronger in the periportal region in females compared to males. Early (1-day) TCPOBOP transcriptional responses were enriched for CAR-bound primary response genes, and for lipogenesis and xenobiotic metabolism and oxidative stress protection pathways; late (2-wk) TCPOBOP responses included many CAR binding-independent secondary response genes, with enrichment for macrophage activation, immune response and cytokine and reactive oxygen species production. Late upstream regulators specific to TCPOBOP-exposed male liver were linked to pro-inflammatory responses and hepatocellular carcinoma progression. TCPOBOP administered weekly to male mice using a high corn oil vehicle activated carbohydrate-responsive transcription factor (MLXIPL)-regulated target genes, dysregulated mitochondrial respiratory and translation regulatory pathways, and induced more advanced liver pathology. Overall, TCPOBOP exposure recapitulates histological and gene expression changes characteristic of emerging steatotic liver disease, including secondary gene responses in liver non-parenchymal cells indicative of transition to a more advanced disease state. Upstream regulators of both the early and late TCPOBOP response genes include novel biomarkers for foreign chemical-induced metabolic dysfunction-associated steatotic liver disease.

The pseudoentropy of allele frequency trajectories, the persistence of variation, and the effective population size.

To concisely describe how genetic variation, at individual loci or across whole genomes, changes over time, and to follow transitory allelic changes, we introduce a quantity related to entropy, that we term pseudoentropy. This quantity emerges in a diffusion analysis of the mean time a mutation segregates in a population. For a neutral locus with an arbitrary number of alleles, the mean time of segregation is generally proportional to the pseudoentropy of initial allele frequencies. After the initial time point, pseudoentropy generally decreases, but other behaviours are possible, depending on the genetic diversity and selective forces present. For a biallelic locus, pseudoentropy and entropy coincide, but they are distinct quantities with more than two alleles. Thus for populations with multiple biallelic loci, the language of entropy suffices. Then entropy, combined across loci, serves as a concise description of genetic variation. We used individual based simulations to explore how this entropy behaves under different evolutionary scenarios. In agreement with predictions, the entropy associated with unlinked neutral loci decreases over time. However, deviations from free recombination and neutrality have clear and informative effects on the entropy's behaviour over time. Analysis of publicly available data of a natural D. melanogaster population, that had been sampled over seven years, using a sliding-window approach, yielded considerable variation in entropy trajectories of different genomic regions. These mostly follow a pattern that suggests a substantial effective population size and a limited effect of positive selection on genome-wide diversity over short time scales.

Plasma growth hormone pulses induce male-biased pulsatile chromatin opening and epigenetic regulation in adult mouse liver.

Sex differences in plasma growth hormone (GH) profiles, pulsatile in males and persistent in females, regulate sex differences in hepatic STAT5 activation linked to sex differences in gene expression and liver disease susceptibility, but little is understood about the fundamental underlying, GH pattern-dependent regulatory mechanisms. Here, DNase-I hypersensitivity site (DHS) analysis of liver chromatin accessibility in a cohort of 18 individual male mice established that the endogenous male rhythm of plasma GH pulse-stimulated liver STAT5 activation induces dynamic, repeated cycles of chromatin opening and closing at several thousand liver DHS and comprises a novel mechanism conferring male bias to liver chromatin accessibility. Strikingly, a single physiological replacement dose of GH given to hypophysectomized male mice restored, within 30 min, liver STAT5 activity and chromatin accessibility at 83% of the dynamic, pituitary hormone-dependent male-biased DHS. Sex-dependent transcription factor binding patterns and chromatin state analysis identified key genomic and epigenetic features distinguishing this dynamic, STAT5-driven mechanism of male-biased chromatin opening from a second GH-dependent mechanism operative at static male-biased DHS, which are constitutively open in male liver. Dynamic but not static male-biased DHS adopt a bivalent-like epigenetic state in female liver, as do static female-biased DHS in male liver, albeit using distinct repressive histone marks in each sex, namely, H3K9me3 at male-biased DHS in female liver and H3K27me3 at female-biased DHS in male liver. Moreover, sex-biased H3K36me3 marks are uniquely enriched at static sex-biased DHS, which may serve to keep these sex-dependent hepatocyte enhancers free of H3K27me3 repressive marks and thus constitutively open. Pulsatile chromatin opening stimulated by endogenous, physiological hormone pulses is thus one of two distinct GH-determined mechanisms for establishing widespread sex differences in hepatic chromatin accessibility and epigenetic regulation, both closely linked to sex-biased gene transcription and the sexual dimorphism of liver function.

Plasma Growth Hormone Pulses Induce Male-biased Pulsatile Chromatin Opening and Epigenetic Regulation in Adult Mouse Liver.

Sex-differences in plasma growth hormone (GH) profiles, pulsatile in males and persistent in females, regulate sex differences in hepatic STAT5 activation linked to sex differences in gene expression and liver disease susceptibility, but little is understood about the fundamental underlying, GH pattern-dependent regulatory mechanisms. Here, DNase hypersensitivity site (DHS) analysis of liver chromatin accessibility in a cohort of 18 individual male mice established that the endogenous male rhythm of plasma GH pulse-stimulated liver STAT5 activation induces dynamic, repeated cycles of chromatin opening and closing at several thousand liver DHS and comprises a novel mechanism conferring male bias to liver chromatin accessibility. Strikingly, a single physiological replacement dose of GH given to hypophysectomized male mice restored, within 30 min, liver STAT5 activity and chromatin accessibility at 83% of the pituitary hormone-dependent dynamic male-biased DHS. Sex-dependent transcription factor binding patterns and chromatin state analysis identified key genomic and epigenetic features distinguishing this dynamic, STAT5-driven mechanism of male-biased chromatin opening from a second GH-dependent mechanism operative at static male-biased DHS, which are constitutively open in male liver. Dynamic but not static male-biased DHS adopt a bivalent-like epigenetic state in female liver, as do static female-biased DHS in male liver, albeit using distinct repressive histone marks in each sex, namely, H3K27me3 at female-biased DHS in male liver, and H3K9me3 at male-biased DHS in female liver. Moreover, sex-biased H3K36me3 marks are uniquely enriched at static sex-biased DHS, which may serve to keep these sex-dependent hepatocyte enhancers free of H3K27me3 repressive marks and thus constitutively open. Pulsatile chromatin opening stimulated by endogenous, physiological hormone pulses is thus one of two distinct GH-determined mechanisms for establishing widespread sex differences in hepatic chromatin accessibility and epigenetic regulation, both closely linked to sex-biased gene transcription and the sexual dimorphism of liver function.

Growth Hormone Accelerates Recovery From Acetaminophen-Induced Murine Liver Injury.

Background and Aims: Acetaminophen (APAP) overdose is the leading cause of acute liver failure, with one available treatment, N-acetyl cysteine (NAC). Yet, NAC effectiveness diminishes about ten hours after APAP overdose, urging for therapeutic alternatives. This study addresses this need by deciphering a mechanism of sexual dimorphism in APAP-induced liver injury, and leveraging it to accelerate liver recovery via growth hormone (GH) treatment. GH secretory patterns, pulsatile in males and near-continuous in females, determine the sex bias in many liver metabolic functions. Here, we aim to establish GH as a novel therapy to treat APAP hepatotoxicity. Approach and Results: Our results demonstrate sex-dependent APAP toxicity, with females showing reduced liver cell death and faster recovery than males. Single-cell RNA sequencing analyses reveal that female hepatocytes have significantly greater levels of GH receptor expression and GH pathway activation compared to males. In harnessing this female-specific advantage, we demonstrate that a single injection of recombinant human GH protein accelerates liver recovery, promotes survival in males following sub-lethal dose of APAP, and is superior to standard-of-care NAC. Alternatively, slow-release delivery of human GH via the safe nonintegrative lipid nanoparticle-encapsulated nucleoside-modified mRNA (mRNA-LNP), a technology validated by widely used COVID-19 vaccines, rescues males from APAP-induced death that otherwise occurred in control mRNA-LNP-treated mice. Conclusions: Our study demonstrates a sexually dimorphic liver repair advantage in females following APAP overdose, leveraged by establishing GH as an alternative treatment, delivered either as recombinant protein or mRNA-LNP, to potentially prevent liver failure and liver transplant in APAP-overdosed patients.

TCDD dysregulation of lncRNA expression, liver zonation and intercellular communication across the liver lobule.

The persistent environmental aryl hydrocarbon receptor agonist and hepatotoxin TCDD (2,3,7,8-tetrachlorodibenzo-p-dioxin) induces hepatic lipid accumulation (steatosis), inflammation (steatohepatitis) and fibrosis. Thousands of liver-expressed, nuclear-localized lncRNAs with regulatory potential have been identified; however, their roles in TCDD-induced hepatoxicity and liver disease are unknown. We analyzed single nucleus (sn)RNA-seq data from control and subchronic (4 wk) TCDD-exposed mouse liver to determine liver cell-type specificity, zonation and differential expression profiles for thousands of lncRNAs. TCDD dysregulated >4000 of these lncRNAs in one or more liver cell types, including 684 lncRNAs specifically dysregulated in liver non-parenchymal cells. Trajectory inference analysis revealed major disruption by TCDD of hepatocyte zonation, affecting >800 genes, including 121 lncRNAs, with strong enrichment for lipid metabolism genes. TCDD also dysregulated expression of >200 transcription factors, including 19 Nuclear Receptors, most notably in hepatocytes and Kupffer cells. TCDD-induced changes in cell-cell communication patterns included marked decreases in EGF signaling from hepatocytes to non-parenchymal cells and increases in extracellular matrix-receptor interactions central to liver fibrosis. Gene regulatory networks constructed from the snRNA-seq data identified TCDD-exposed liver network-essential lncRNA regulators linked to functions such as fatty acid metabolic process, peroxisome and xenobiotic metabolism. Networks were validated by the striking enrichments that predicted regulatory lncRNAs showed for specific biological pathways. These findings highlight the power of snRNA-seq to discover functional roles for many xenobiotic-responsive lncRNAs in both hepatocytes and liver non-parenchymal cells and to elucidate novel aspects of foreign chemical-induced hepatotoxicity and liver disease, including dysregulation of intercellular communication within the liver lobule.

Direct and systemic actions of growth hormone receptor (GHR)-signaling on hepatic glycolysis, de novo lipogenesis and insulin sensitivity, associated with steatosis.

BACKGROUND: Evidence is accumulating that growth hormone (GH) protects against the development of steatosis and progression of non-alcoholic fatty liver disease (NAFLD). GH may control steatosis indirectly by altering systemic insulin sensitivity and substrate delivery to the liver and/or by the direct actions of GH on hepatocyte function. APPROACH: To better define the hepatocyte-specific role of GH receptor (GHR) signaling on regulating steatosis, we used a mouse model with adult-onset, hepatocyte-specific GHR knockdown (aHepGHRkd). To prevent the reduction in circulating insulin-like growth factor 1 (IGF1) and the subsequent increase in GH observed after aHepGHRkd, subsets of aHepGHRkd mice were treated with adeno-associated viral vectors (AAV) driving hepatocyte-specific expression of IGF1 or a constitutively active form of STAT5b (STAT5bCA). The impact of hepatocyte-specific modulation of GHR, IGF1 and STAT5b on carbohydrate and lipid metabolism was studied across multiple nutritional states and in the context of hyperinsulinemic:euglycemic clamps. RESULTS: Chow-fed male aHepGHRkd mice developed steatosis associated with an increase in hepatic glucokinase (GCK) and ketohexokinase (KHK) expression and de novo lipogenesis (DNL) rate, in the post-absorptive state and in response to refeeding after an overnight fast. The aHepGHRkd-associated increase in hepatic KHK, but not GCK and steatosis, was dependent on hepatocyte expression of carbohydrate response element binding protein (ChREBP), in re-fed mice. Interestingly, under clamp conditions, aHepGHRkd also increased the rate of DNL and expression of GCK and KHK, but impaired insulin-mediated suppression of hepatic glucose production, without altering plasma NEFA levels. These effects were normalized with AAV-mediated hepatocyte expression of IGF1 or STAT5bCA. Comparison of the impact of AAV-mediated hepatocyte IGF1 versus STAT5bCA in aHepGHRkd mice across multiple nutritional states, indicated the restorative actions of IGF1 are indirect, by improving systemic insulin sensitivity, independent of changes in the liver transcriptome. In contrast, the actions of STAT5b are due to the combined effects of raising IGF1 and direct alterations in the hepatocyte gene program that may involve suppression of BCL6 and FOXO1 activity. However, the direct and IGF1-dependent actions of STAT5b cannot fully account for enhanced GCK activity and lipogenic gene expression observed after aHepGHRkd, suggesting other GHR-mediated signals are involved. CONCLUSION: These studies demonstrate hepatocyte GHR-signaling controls hepatic glycolysis, DNL, steatosis and hepatic insulin sensitivity indirectly (via IGF1) and directly (via STAT5b). The relative contribution of these indirect and direct actions of GH on hepatocytes is modified by insulin and nutrient availability. These results improve our understanding of the physiologic actions of GH on regulating adult metabolism to protect against NAFLD progression.

Dysregulation of murine long noncoding single-cell transcriptome in nonalcoholic steatohepatitis and liver fibrosis.

LncRNAs comprise a heterogeneous class of RNA-encoding genes typified by low expression, nuclear enrichment, high tissue-specificity, and functional diversity, but the vast majority remain uncharacterized. Here, we assembled the mouse liver noncoding transcriptome from >2000 bulk RNA-seq samples and discovered 48,261 liver-expressed lncRNAs, a majority novel. Using these lncRNAs as a single-cell transcriptomic reference set, we elucidated lncRNA dysregulation in mouse models of high fat diet-induced nonalcoholic steatohepatitis and carbon tetrachloride-induced liver fibrosis. Trajectory inference analysis revealed lncRNA zonation patterns across the liver lobule in each major liver cell population. Perturbations in lncRNA expression and zonation were common in several disease-associated liver cell types, including nonalcoholic steatohepatitis-associated macrophages, a hallmark of fatty liver disease progression, and collagen-producing myofibroblasts, a central feature of liver fibrosis. Single-cell-based gene regulatory network analysis using bigSCale2 linked individual lncRNAs to specific biological pathways, and network-essential regulatory lncRNAs with disease-associated functions were identified by their high network centrality metrics. For a subset of these lncRNAs, promoter sequences of the network-defined lncRNA target genes were significantly enriched for lncRNA triplex formation, providing independent mechanistic support for the lncRNA-target gene linkages predicted by the gene regulatory networks. These findings elucidate liver lncRNA cell-type specificities, spatial zonation patterns, associated regulatory networks, and temporal patterns of dysregulation during hepatic disease progression. A subset of the liver disease-associated regulatory lncRNAs identified have human orthologs and are promising candidates for biomarkers and therapeutic targets.

TCDD dysregulation of lncRNA expression, liver zonation and intercellular communication across the liver lobule.

The persistent environmental aryl hydrocarbon receptor agonist and hepatotoxin TCDD (2,3,7,8-tetrachlorodibenzo- p -dioxin) induces hepatic lipid accumulation (steatosis), inflammation (steatohepatitis) and fibrosis. Thousands of liver-expressed, nuclear-localized lncRNAs with regulatory potential have been identified; however, their roles in TCDD-induced hepatoxicity and liver disease are unknown. We analyzed single nucleus (sn)RNA-seq data from control and chronic TCDD-exposed mouse liver to determine liver cell-type specificity, zonation and differential expression profiles for thousands of IncRNAs. TCDD dysregulated >4,000 of these lncRNAs in one or more liver cell types, including 684 lncRNAs specifically dysregulated in liver non-parenchymal cells. Trajectory inference analysis revealed major disruption by TCDD of hepatocyte zonation, affecting >800 genes, including 121 IncRNAs, with strong enrichment for lipid metabolism genes. TCDD also dysregulated expression of >200 transcription factors, including 19 Nuclear Receptors, most notably in hepatocytes and Kupffer cells. TCDD-induced changes in cellâ€"cell communication patterns included marked decreases in EGF signaling from hepatocytes to non-parenchymal cells and increases in extracellular matrix-receptor interactions central to liver fibrosis. Gene regulatory networks constructed from the snRNA-seq data identified TCDD-exposed liver network-essential lncRNA regulators linked to functions such as fatty acid metabolic process, peroxisome and xenobiotic metabolic. Networks were validated by the striking enrichments that predicted regulatory IncRNAs showed for specific biological pathways. These findings highlight the power of snRNA-seq to discover functional roles for many xenobiotic-responsive lncRNAs in both hepatocytes and liver non-parenchymal cells and to elucidate novel aspects of foreign chemical-induced hepatotoxicity and liver disease, including dysregulation of intercellular communication within the liver lobule.

Sex matters in liver fat regulation.

Growth hormone orchestrates a complex, sex-dependent balancing act.

Type-I interferon signaling is essential for robust metronomic chemo-immunogenic tumor regression in murine breast cancer.

Many patients with breast cancer have a poor prognosis with limited therapeutic options. Here, we investigated the potential of chemo-immunogenic therapy as an avenue of treatment. We utilized two syngeneic mouse mammary tumor models, 4T1 and E0771, to examine the chemo-immunogenic potential of cyclophosphamide and the mechanistic contributions of cyclophosphamide-activated type-I interferon (IFN) signaling to therapeutic activity. Chemically-activated cyclophosphamide induced robust IFNα/ß receptor-1-dependent signaling linked to hundreds of IFN-stimulated gene responses in both cell lines. Further, in 4T1 tumors, cyclophosphamide given on a medium-dose, 6-day intermittent metronomic schedule induced strong IFN signaling but comparatively weak immune cell infiltration associated with long-term tumor growth stasis. Induction of IFN signaling was somewhat weaker in E0771 tumors but was followed by widespread downstream gene responses, robust immune cell infiltration and extensive, prolonged tumor regression. The immune dependence of these effective anti-tumor responses was established by CD8 T-cell immunodepletion, which blocked cyclophosphamide-induced E0771 tumor regression and led to tumor stasis followed by regrowth. Strikingly, IFNα/ß receptor-1 antibody blockade was even more effective in preventing E0771 immune cell infiltration and blocked the major tumor regression induced by cyclophosphamide treatment. Type-I IFN signaling is thus essential for the robust chemo-immunogenic response of these tumors to cyclophosphamide administered on a metronomic schedule.

Long non-coding RNA G23Rik attenuates fasting-induced lipid accumulation in mouse liver.

Peroxisome proliferator-activated receptor α (PPARα) is a key mediator of lipid metabolism and metabolic stress in the liver. A recent study revealed that PPARα-dependent long non-coding RNAs (lncRNAs) play an important role in modulating metabolic stress and inflammation in the livers of fasted mice. Here hepatic lncRNA 3930402G23Rik (G23Rik) was found to have active peroxisome proliferator response elements (PPREs) within its promoter and is directly regulated by PPARα. Although G23Rik RNA was expressed to varying degrees in several tissues, the PPARα-dependent regulation of this lncRNA was only observed in the liver. Pharmacological activation of PPARα induced PPARα recruitment at the G23Rik promoter and a pronounced increase in hepatic G23Rik lncRNA expression. A G23Rik-null mouse line was developed to further characterize the function of this lncRNA in the liver. G23Rik-null mice were more susceptible to hepatic lipid accumulation in response to acute fasting. Histological analysis further revealed a pronounced buildup of lipid droplets and a significant increase in neutral triglycerides and lipids as indicated by enhanced oil red O staining of liver sections. Hepatic cholesterol, non-esterified fatty acid, and triglyceride levels were significantly elevated in G23Rik-null mice and associated with induction of the lipid-metabolism related gene Cd36. These findings provide evidence for a lncRNA dependent mechanism by which PPARα attenuates hepatic lipid accumulation in response to metabolic stress through lncRNA G23Rik induction.

Interplay Between GH-regulated, Sex-biased Liver Transcriptome and Hepatic Zonation Revealed by Single-Nucleus RNA Sequencing.

The zonation of liver metabolic processes is well-characterized; however, little is known about the cell type-specificity and zonation of sexually dimorphic gene expression or its growth hormone (GH)-dependent transcriptional regulators. We address these issues using single-nucleus RNA-sequencing of 32 000 nuclei representing 9 major liver cell types. Nuclei were extracted from livers from adult male and female mice; from males infused with GH continuously, mimicking the female plasma GH pattern; and from mice exposed to TCPOBOP, a xenobiotic agonist ligand of the nuclear receptor CAR that perturbs sex-biased gene expression. Analysis of these rich transcriptomic datasets revealed the following: 1) expression of sex-biased genes and their GH-dependent transcriptional regulators is primarily restricted to hepatocytes and is not a feature of liver nonparenchymal cells; 2) many sex-biased transcripts show sex-dependent zonation within the liver lobule; 3) gene expression is substantially feminized both in periportal and pericentral hepatocytes when male mice are infused with GH continuously; 4) sequencing nuclei increases the sensitivity for detecting thousands of nuclear-enriched long-noncoding RNAs (lncRNAs) and enables determination of their liver cell type-specificity, sex-bias and hepatocyte zonation profiles; 5) the periportal to pericentral hepatocyte cell ratio is significantly higher in male than female liver; and 6) TCPOBOP exposure disrupts both sex-specific gene expression and hepatocyte zonation within the liver lobule. These findings highlight the complex interconnections between hepatic sexual dimorphism and zonation at the single-cell level and reveal how endogenous hormones and foreign chemical exposure can alter these interactions across the liver lobule with large effects both on protein-coding genes and lncRNAs.

Theoretically quantifying the direct and indirect benefits of vaccination against SARS-CoV-2 in terms of avoided deaths.

The Coronavirus Disease 2019 (COVID)-19 pandemic has placed unprecedented pressures on societies around the world. Successful vaccines, developed against the spike protein of the Severe Acute Respiratory Syndrome Virus 2 (SARS-CoV-2) virus, offer hope that new hospitalisations and new deaths will subside. However, vaccination takes place in a dynamic environment. For example, new variants of the disease may occur where the effectiveness of a vaccine lies below that of the original target of the vaccine, while changes in the behaviour of a population are accompanied by a changed basic reproduction number. Here, we aim to understand how changes in values of basic parameters affect the benefits of vaccination at the direct level, of the individuals vaccinated, and at the indirect level, of the wider, unvaccinated community. We work within the framework of a Susceptible-Infected-Recovered model, and produce a metric for the benefits of vaccination, at both direct and indirect levels, in terms of the number of avoided deaths. Taking into account the initial prevalence of a SARS-CoV-2 infection, the mortality rate of the disease, the basic reproduction number, the vaccination rate, and the effectiveness of a vaccine, we explore how these basic parameters affect the benefits of vaccination. We find a range of situations where indirect benefits of vaccination outweigh direct benefits. This especially occurs at lower rates of vaccination (20% - [Formula: see text]) and intermediate values of the basic reproduction number (1-1.5). The indirect benefits can be substantial, in some cases being more than 400% of the direct benefits. For an initial prevalence of SARS-CoV-2 infection of 2%, a basic reproduction number of 1.2, a mortality rate of 2%, and a vaccine effectiveness of 95%, our findings show, for a population of 500,000 people, where 100,000 susceptible individuals are vaccinated, that approximately 2200 deaths are avoided. However, approximately 600 of these deaths are avoided amongst vaccinated individuals, while approximately 1600 deaths are avoided in the wider, unvaccinated community.

Spatial frequency domain imaging for monitoring immune-mediated chemotherapy treatment response and resistance in a murine breast cancer model.

Spatial Frequency Domain Imaging (SFDI) can provide longitudinal, label-free, and widefield hemodynamic and scattering measurements of murine tumors in vivo. Our previous work has shown that the reduced scattering coefficient (µ's) at 800 nm, as well as the wavelength dependence of scattering, both have prognostic value in tracking apoptosis and proliferation during treatment with anti-cancer therapies. However, there is limited work in validating these optical biomarkers in clinically relevant tumor models that manifest specific treatment resistance mechanisms that mimic the clinical setting. It was recently demonstrated that metronomic dosing of cyclophosphamide induces a strong anti-tumor immune response and tumor volume reduction in the E0771 murine breast cancer model. This immune activation mechanism can be blocked with an IFNAR-1 antibody, leading to treatment resistance. Here we present a longitudinal study utilizing SFDI to monitor this paired responsive-resistant model for up to 30 days of drug treatment. Mice receiving the immune modulatory metronomic cyclophosphamide schedule had a significant increase in tumor optical scattering compared to mice receiving cyclophosphamide in combination with the IFNAR-1 antibody (9% increase vs 10% decrease on day 5 of treatment, p < 0.001). The magnitude of these differences increased throughout the duration of treatment. Additionally, scattering changes on day 4 of treatment could discriminate responsive versus resistant tumors with an accuracy of 78%, while tumor volume had an accuracy of only 52%. These results validate optical scattering as a promising prognostic biomarker that can discriminate between treatment responsive and resistant tumor models.

Constitutively Active STAT5b Feminizes Mouse Liver Gene Expression.

STAT5 is an essential transcriptional regulator of the sex-biased actions of GH in the liver. Delivery of constitutively active STAT5 (STAT5CA) to male mouse liver using an engineered adeno-associated virus with high tropism for the liver is shown to induce widespread feminization of the liver, with extensive induction of female-biased genes and repression of male-biased genes, largely mimicking results obtained when male mice are given GH as a continuous infusion. Many of the STAT5CA-responding genes were associated with nearby (< 50 kb) sites of STAT5 binding to liver chromatin, supporting the proposed direct role of persistently active STAT5 in continuous GH-induced liver feminization. The feminizing effects of STAT5CA were dose-dependent; moreover, at higher levels, STAT5CA overexpression resulted in some histopathology, including hepatocyte hyperplasia, and increased karyomegaly and multinuclear hepatocytes. These findings establish that the persistent activation of STAT5 by GH that characterizes female liver is by itself sufficient to account for the sex-dependent expression of a majority of hepatic sex-biased genes. Moreover, histological changes seen when STAT5CA is overexpressed highlight the importance of carefully evaluating such effects before considering STAT5 derivatives for therapeutic use in treating liver disease.

Impact of Neonatal Activation of Nuclear Receptor CAR (Nr1i3) on Cyp2 Gene Expression in Adult Mouse Liver.

Perinatal exposure to environmental chemicals is proposed to reprogram development and alter disease susceptibility later in life. Supporting this, neonatal activation of the nuclear receptor constitutive androstane receptor (CAR) (Nr1i3) by TCPOBOP was previously reported to induce persistent expression of mouse hepatic Cyp2 genes into adulthood, and was attributed to long-term epigenetic memory of the early life exposure. Here, we confirm that the same high-dose neonatal TCPOBOP exposure studied previously (3 mg/kg, 15x ED50) does indeed induce prolonged (12 weeks) increases in hepatic Cyp2 expression; however, we show that the persistence of expression can be fully explained by the persistence of residual TCPOBOP in liver tissue. When the long-term presence of TCPOBOP in tissue was eliminated by decreasing the neonatal TCPOBOP dose 22-fold (0.67× ED50), strong neonatal increases in hepatic Cyp2 expression were still obtained but did not persist into adulthood. Furthermore, the neonatal ED50-range TCPOBOP exposure did not sensitize mice to a subsequent, low-dose TCPOBOP treatment. In contrast, neonatal treatment with phenobarbital, a short half-life (t1/2 = 8 h) agonist of CAR and PXR (Nr1i2), induced high-level neonatal activation of Cyp2 genes and also altered their responsiveness to low-dose phenobarbital exposure at adulthood by either increasing (Cyp2b10) or decreasing (Cyp2c55) expression. Thus, neonatal xenobiotic exposure can reprogram hepatic Cyp2 genes and alter their responsiveness to exposures later in life. These findings highlight the need to carefully consider xenobiotic dose, half-life, and persistence in tissue when evaluating the long-term effects of early life environmental chemical exposures.

Correcting Bias in Allele Frequency Estimates Due to an Observation Threshold: A Markov Chain Analysis.

There are many problems in biology and related disciplines involving stochasticity, where a signal can only be detected when it lies above a threshold level, while signals lying below threshold are simply not detected. A consequence is that the detected signal is conditioned to lie above threshold, and is not representative of the actual signal. In this work, we present some general results for the conditioning that occurs due to the existence of such an observational threshold. We show that this conditioning is relevant, for example, to gene-frequency trajectories, where many loci in the genome are simultaneously measured in a given generation. Such a threshold can lead to severe biases of allele frequency estimates under purifying selection. In the analysis presented, within the context of Markov chains such as the Wright-Fisher model, we address two key questions: (1) "What is a natural measure of the strength of the conditioning associated with an observation threshold?" (2) "What is a principled way to correct for the effects of the conditioning?". We answer the first question in terms of a proportion. Starting with a large number of trajectories, the relevant quantity is the proportion of these trajectories that are above threshold at a later time and hence are detected. The smaller the value of this proportion, the stronger the effects of conditioning. We provide an approximate analytical answer to the second question, that corrects the bias produced by an observation threshold, and performs to reasonable accuracy in the Wright-Fisher model for biologically plausible parameter values.

Increased brain volume from higher cereal and lower coffee intake: shared genetic determinants and impacts on cognition and metabolism.

It is unclear how different diets may affect human brain development and if genetic and environmental factors play a part. We investigated diet effects in the UK Biobank data from 18,879 healthy adults and discovered anticorrelated brain-wide gray matter volume (GMV)-association patterns between coffee and cereal intake, coincidence with their anticorrelated genetic constructs. The Mendelian randomization approach further indicated a causal effect of higher coffee intake on reduced total GMV, which is likely through regulating the expression of genes responsible for synaptic development in the brain. The identified genetic factors may further affect people's lifestyle habits and body/blood fat levels through the mediation of cereal/coffee intake, and the brain-wide expression pattern of gene CPLX3, a dedicated marker of subplate neurons that regulate cortical development and plasticity, may underlie the shared GMV-association patterns among the coffee/cereal intake and cognitive functions. All the main findings were successfully replicated. Our findings thus revealed that high-cereal and low-coffee diets shared similar brain and genetic constructs, leading to long-term beneficial associations regarding cognitive, body mass index (BMI), and other metabolic measures. This study has important implications for public health, especially during the pandemic, given the poorer outcomes of COVID-19 patients with greater BMIs.

Harnessing natural variation to identify cis regulators of sex-biased gene expression in a multi-strain mouse liver model.

Sex differences in gene expression are widespread in the liver, where many autosomal factors act in tandem with growth hormone signaling to regulate individual variability of sex differences in liver metabolism and disease. Here, we compare hepatic transcriptomic and epigenetic profiles of mouse strains C57BL/6J and CAST/EiJ, representing two subspecies separated by 0.5-1 million years of evolution, to elucidate the actions of genetic factors regulating liver sex differences. We identify 144 protein coding genes and 78 lncRNAs showing strain-conserved sex bias; many have gene ontologies relevant to liver function, are more highly liver-specific and show greater sex bias, and are more proximally regulated than genes whose sex bias is strain-dependent. The strain-conserved genes include key growth hormone-dependent transcriptional regulators of liver sex bias; however, three other transcription factors, Trim24, Tox, and Zfp809, lose their sex-biased expression in CAST/EiJ mouse liver. To elucidate the observed strain specificities in expression, we characterized the strain-dependence of sex-biased chromatin opening and enhancer marks at cis regulatory elements (CREs) within expression quantitative trait loci (eQTL) regulating liver sex-biased genes. Strikingly, 208 of 286 eQTLs with strain-specific, sex-differential effects on expression were associated with a complete gain, loss, or reversal of the sex differences in expression between strains. Moreover, 166 of the 286 eQTLs were linked to the strain-dependent gain or loss of localized sex-biased CREs. Remarkably, a subset of these CREs apparently lacked strain-specific genetic variants yet showed coordinated, strain-dependent sex-biased epigenetic regulation. Thus, we directly link hundreds of strain-specific genetic variants to the high variability in CRE activity and expression of sex-biased genes and uncover underlying genetically-determined epigenetic states controlling liver sex bias in genetically diverse mouse populations.